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Recently, a new approach, based on boundary conformal field theory, has been applied to a variety 
of quantum impurity problems in condensed matter and particle physics. A particularly enlightening 
example is the multi-channel Kondo problem. In this review some earlier approaches to the Kondo 
problem are discussed, the needed material on boundary conformal field theory is developed and 
then this new method is applied to the multi-channel Kondo problem. 
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I. RENORMALIZATION GROUP AND FERMI LIQUID APPROACHES TO THE KONDO EFFECT 



A. Introduction to the Kondo Effect 



Most mechanisms contributing to the resistivity of metals, p(T), give either p(T) decreasing to 0, 
as T — > (phonons or electron-electron interactions), or p(T) — > constant, as T — ► (non-magnetric 
impurities). However, metals ppntaining magnetic impurities show a p(T) which increases as T — > 0. 
This was explained by KondoEJ in 1964 using a simple Hamiltonian: 

H = + xs- 5>jf v-p (i-i) 

kct kk' 

where i>% a 's are conduction electron annihilation operators, (of momentum k, spin a) and S repre- 
sents the spin of the magnetic impurity with 

[S a ,S b ] = ie abc S c . 

The interaction term represents an impurity spin interacting with the electron spin at x = 0. 

With the above Hamiltonian, the Born approximation gives: p(T) ~ A 2 , independent of T. The 
next order term has a divergent coefficient at T = 0: 

p(T) ~ [A + u\ 2 In ^ + ...] 2 (1.2) 

Here D is the band-width, v the density of states. This result stimulated an enormous amount of 
theoretical work. As pNozieres put it, "Theorists 'diverged' on their own, leaving the experiment 
realities way behind" fl What happens at low T, i.e. T ~ Tk = De~~? In that case the 0(A 2 ) 
term will be as big as the term of O(A). What about the 0(A 3 ) term? Such questions helped lead 
to the development of the renormalization group needed to understand the problem. 

In particle physics such a growth of a coupling constant at low energies explains quark-icon- 
finement (1973) and "asymptotic freedom" at E — > oo. To solve these problems, Wilsonu de- 
veloped a very powerful numejeieal renormalization group approach. The Kondo model was aJsi 
"solved" by the Bethe ansatzu'Q which gives the specific heat and magnetization. Nozierespl 
following ideas of Andersonu and Wilson,Q developed a very simple, and in a sense exact, pic- 
ture of th e il ow | C-| bqha y4pur. With A. Ludwig, I have generalized and reformulated Nozieres' 
approacrJMeBElma using recent results in conformal field theory. The latter approach 
is very general and can be annUed-.to,.-a| jmmber of other problems including multj-cbannel 
and higher spin Kondo Bffect,B't^BtM ! Bn two (or more)-impurity Kondo effects, LD'Efl im- 
purity assisted t un it jiol ingjlj impurities in one-dimensional conductors ("quantum wires" )E3 or ID 
antiferromagnets^t 2 ] baryon-monopole interaction,E3. Some of these problems, including the 
multi-channel Kondo effect, exhibit non- Fermi liquid behaviour. These are among the very few 
exactly solved problems that do this (the others are ID Luttinger liquids). It has been sug- 
gested that this paa* be connected with exotic behaviour of certain compounds, including high-T c 



superconductors.l 



B. Renormalization Group 

We could integrate out ip(k) for k far from kp, the Fermi wave-vector, and successively reduce the 
band- width D to obtain a new effective interaction. [See Figure (Q).] This is hard to do exactly. At 
weak coupling one can do it perturbatively in A. With the simplest approach, real-time, time-ordered 
perturbation theory, we expand 



Texp 



i\ / S(t) ■ ip^ip(0,t) 



where the fields are in the interaction picture. 
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FIG. 1. Reduction of the cut-off from D to D'. 

As S(t) is independent of t, we simply multiply powers of 5* using 

[S a ,S b ] = ie abc S", S 2 = s(s + 1). 

We must time-order S's which don't commute. The first few diagrams are shown in Figure (||). In 
2nd order in A, we have: 

_^ J dtdt'T{S a {t)S\t'))-T[^\t)^m^^)^{t% 
which can be reduced, using Wick's theorem, to: 



\ z I dt dt ip 



'„/,t 



2 ' 2 



^T(ip{t)^{t')){0{t - t')S a S b + 0(t' - t)S b S a ) 



X 2 



— / dt dt'tp 1 -ip- S sn(t-t'){^(i)ip \t')) 



2 J "2 

where sn (t — t ) is the sign- function which arises from T-ordering spins. 



(1.3) 
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FIG. 2. Feynman diagrams contributing to renormalization of the Kondo coupling constant to third order. 
We see that the integral 

dt 



dte(t)G(t) = -i J ~ 

is divergent in the infrared limit: t — > oo, where G(i) = (^(t)ip^ (0)}. But we only integrate out 
electrons with D' < k < D which gives In D/D' . To do it explicitly we use the Fourier transformed 
form: 



f d 3 k 


f du> 


r l 


1 WJ 


' 2^ 


,uo + S 



d A k 1 



(2 7 r)3|e fe | - 2 ^ D , 



lu — efc + i<5sn(et 



— = 2^ In — . 



Thus 



and 



<SA = vX 2 In-^-, 



rfA 2 
= — vX . 



d\nD 

We see that lowering the band cut-off increases A or, defining a length-dependent cut-off, I ~ vf/D, 



dlnl 



Integrating the equation (equivalent to performing an infinite sum of diagrams), gives: 

Ac 



1 _ u \ i n 



If Ao > (antiferromagnetic), then A e ff(D) diverges at D ~ Tk ~ Doe vA o , If Ao < (ferromag- 
netic), A eff (D) -> 0. See Figure (§). 



(1.4) 



(1.5) 
(1.6) 

(1.7) 
(1.8) 

(1.9) 

(1.10) 







eff 



FIG. 3. RG flow of the Kondo coupling. 



The behaviour at temperature T is determined by X e g(T): p(T) — > as T — > for the ferromag- 
netic case. What happens for the antiferromagnetic case? 
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C. Mapping to a One-Dimensional Model 

The above discussion can be simplified if we map the model into a one dimensional one. We 
assume a spherically symmetric e(k), 

k 2 

<k) = - e F ^v F {k-k F ), (1.11) 

Am 

and a 5— function Kondo interaction. There is only s-wave scattering, i.e. 

ip(k) = tpo (fc) + higher harmonics, 

V47rfc 



Ho = J dkipQ k ipok(-{k) + higher harmonics, 



a 



#INT = XvpV J dkdk ^o,k2^°- k ' ' S ' ( l - 12 ) 

where v = k F /2ir 2 VF is the density of states per spin. This can also be written in terms of radial 
co-ordinate. We eliminate all modes except for a band width 2D: \k — k F \ < D. Defining left and 
right movers (incoming and outgoing waves), 

<f L , R (r)= f dke ±lkr Mk + k F ), => ^(0)=Vfl(0), (1.13) 

J —A 

we have 

H — — / drOipti-^-ipL — i'ni-T-ipR) (note the unconventional normalization) , 
Z7r J dr dr 

H mT = v F \^ L (Q)^4, L (0)-S. (1.14) 

Here we have redefined a dimensionless Kondo coupling, A — > \v. Using the notation 

tpL = tPl(x,t) = 4>l(z = r + ix), i/)r(x,t) = ip R (z* = r — ix), (1-15) 

where r is imaginary time and x = r, (and we set v F = 1) we have 

(M*)1>t(o)) = \, <iM**)i&(o)> = y*- (L16) 

Alternatively, since 

V>z,(0,t) = Vh(0,t) ^ L =i> L {,z), iPr = ^ r {z*), (1.17) 

we may consider tpR to be the continuation of ipL to the negative r-axis: 

ip R (x,T) = ^ l (-x,t). (1-18) 

Now we obtain a relativistic (1 + 1) dimensional field theory ( a "crural" one, containing left-movers 
only) interacting with the impurity at x = with 



VF f°° , ,t . d 
2tt ( 

and #int as in Eq. ( |L14| ), See Figure (g) 



#o = / dxrl)li—il)L (1-19) 
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FIG. 4. Reflecting the left-movers to the negative axis. 

D. Fermi Liquid Approach at Low T 

What is the T — > behavior of the antiferromagetic Kondo model? The simplest assumption is 
A c jj — > oo. But what does that really mean? Consider the strong coupling limit of a lattice model,0 
for convenience, in spatial dimension D = 1. (D doesn't really matter since we can always reduce 
the model to D — 1.) 

H = t ^(VjVi+i + 1>l +1 lH) + ■ V-Jf^o (1.20) 

i 

Consider the limit A >> tj. The groundstate of the interaction term will be the following con- 
figuration: one electron at the site forms a singlet with the impurity: | ~{[[) — | JJ-f). (We as- 
sume Simp = 1/2). Now we do perturbation theory in t. We have the following low energy states: 
an arbitary electron configuration occurs on all other sites-but other electrons or holes are forbidden 
to enter the site-0, since that would destroy the singlet state, costing an energy, AE ~ A >> t. 
Thus we simply form free electron Bloch states with the boundary condition (f)(0) = 0, where cp(i) is 
the single-electron wave-function. Note that at zero Kondo coupling, the parity even single particle 
wave-functions are of the form <f>(i) = cos ki and the parity odd ones are of the form (j>(i) = sin ki. 
On the other hand, at A — *• oo the parity even wave- functions become <j>(i) = |sinfci|, while the 
parity odd ones are unaffected. 

The behaviour of the parity even channel corresponds to a n/2 phase shift in the s-wave channel. 

~ e ~ im +e + ™e im , S = tt/2. (1.21) 
In terms of left and right movers on r > we have changed the boundary condition, 

i>L(o) = i> R (o), a = o, 

i> L (0) = -i> R (Q), A = oo. (1.22) 

The strong coupling fixed point is the same as the weak coupling fixed point except for a change 
in boundary conditions (and the removal of the impurity). In terms of the left-moving description 
of the P-even sector, the phase of the left-mover is shifted by n as it passes the origin. Imposing 
another boundary condition a distance I away quantizes k: 

il>{i) = Vi(0 + tMO = <M0 + i>L(-i) = o, 

A = 0: fe=~(n+l/2) 

A = oo: fc=^ (1.23) 



G 



Near the Fermi surface the energies are linearly spaced. Assuming particle-hole symmetry, the 
Fermi energy lies midway between levels or on a level. [See Figures (^) and (|^).] The two situations 
switch with the phase shift. Wilson's numerical RG schemeu involves calculating the low- lying 
spectrum numerically and looking for this shift. This indicates that A renormalizes to oo even if it 
is initially small. However, now we expect the screening to take place over a longer length scale 



V£ VF 1/„A 

T k ~ D B ■ 



In other words, the wave function of the screening electron has this scale. We get low energy Bloch 
states of free electrons only for fc — kp\ << l/£ (so we must take I » £). [See Figure (R).] The free 
electron theory with a phase shift corresponds to a universal stable low energy fixed point for the 
Kondo problem. This observation determines the T — resistivity for an array of Kondo impurities 
at random locations of low density m. It is the same as for non-magnetic s-wave scatterers with a 
7r/2 phase shift at the Fermi energy. S — tv/2 gives the so-called unitary limit resistivity: 



Pu = 



3m 



(1.24) 



(1.25) 




k-k 



FIG. 5. Free fermion energy levels with antiperiodic boundary conditions. 
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FIG. 6. Free fermion energy levels with periodic boundary conditions. 

screening cloud Bloch wave-function 




J_ 
k-k, 

FIG. 7. Non-interacting Bloch states with a vanishing boundary condition occur for \k — fcjr| << vf/Tk- 

The low-T behaviour, so far, seems trivial. Much of the interesting behaviour comes from the 
leading irrelevant operator. The impurity spin has disappeared (screened) from the description of 
the low -T physics. However certain interactions between electrons are generated (at the impurity 
site only) in the process of eliminating the impurity spin. We can determine these by simply writing 
the lowest dimension operators allowed by symmetry. 

It is simplest to work in the ID formulation, with left-movers only. We write the interaction in 
terms of ipL, obeying the new boundary condition (but not the impurity spin). The dimension of 
the operator is determined as in ID field theory 

H = J dx^ L i^L + .... (1.26) 

The length and time dimensions are equivalent (we convert with vf), 

[H] = E=> [ip] = 
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The interactions are local 



MT = 5^0^ = 0), [Ai] + [Oi] = 1. 



So \i has negative energy dimension if [Oi] > 1, implying that it is irrelevant. In RG theory one 
usually defines a dimensionless coupling constant by multiplying powers of the cut-off D, if 

= E~ a , A; ee \ t D a , 

Ai decreases as we lower D: 

Such a coupling, with a > 0, produces no infrared divergences in perturbation theory. The ultraviolet 
ones are cancelled by the expicit factors of the ultraviolet cut-off, D, appearing in the Lagrangian, 
XO/D a . What are the lowest dimension operators allowed by symmetry? Consider ip +a (0)ip a (0). 
This has d=l. However, it is not allowed because it breaks particle-hole symmetry. If particle-hole 
symmetry is broken then we do get this, a potential scattering term; it adds a term to the phase 
shift. Consider another term, 

^"^"i )-*^"^ )- (L28) 

This has d—2 . This term produces a fc-dependent phase shift. The only other term with d < 2 is 
ip^ip^tpl. This term represents the electron-electron interaction induced by an impurity spin- 
flip. The first electron flips the impurity spin. This makes it possible for the second electron to flip 
it back if the electron spin is correct. These are the only d < 2 operators. There are no relevant 
(d < 1) operators, implying the stability of the low energy fixed point. Note that, by contrast, the 
high energy, zero Kondo coupling, fixed point is unstable. The dimension 1 operator (the Kondo 
interaction) can occur there because of the presence of the impurity spin. 

We can't calculate these two coupling constants exactly except by using complicated methods: 
Wilson's numerical melthod or the Bethe ansatz. They both have dimension E~ x . We expect them 
to be 0[1/Tk] by a standard scaling argument. That is, functions of the cut-off D and coupling 
constant A can be replaced by functions of the reduced cut-off, D' and the renormalized coupling 
constant, A en -(D'): f[D, A] = f[D', A en -(D')]. We can lower the cut-off down to Tk where A e g- is 
O(l) so / = /(Tk, 1) = /(Tk). This is a characteristic scale introduced by the infrared divergences 
of perturbation theory. For Tk << D(\ << 1), Nozieres argued that the two irrelevant coupling 
constants have a universal ratio. So there is only one unknown parameter ("the Wilson number"). 
Essentially all low-temperature information is given by this irrelevant coupling constant, if it was 
not already determined by the 7r/2 phase shift at Uf- I will give a different derivation of the ratio 
of the two coupling constants later using conformal field theory. We now simply do perturbation 
theory in the irrelevant coupling constant ~ 1/Tk- We can determine powers of T by dimensional 
analysis. For the specific heat we find: 

c -ir F lT+a h (L29) 

This is the specific heat for the one-dimensional system with a single impurity at the origin. Note 
that the first term is simply the specific heat of the free system, proportional to system length. The 
second term is independent of length and is the impurity specific heat. It is the result of first order 
perturbation theory in the irrelevant coupling constant, of 0(1/Tk)- The (linear) power of T can be 
fixed by dimensional analysis; a is a pure number. Note that while this is formally an "irrelevant" 
contribution, it in fact gives the leading impurity specific heat at low T. To obtain the specific heat 
for the three dimensional system we simply multiply the first term by the ratio vV/(1/2ttvf) i.e. the 
ratio of densities of states per unit energy. For a dilute random array the last term gets multiplied 
by the number of impurities. At high T we get, approximately, the entropy for a decoupled s=l/2 
impurity: 

S(T) = -^-T + ln2. (1.30) 

6Vf 

At low T, the impurity entropy decreases to 0: 
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Av F Ik 



(1.31) 



In general we may write: 

S(T)-^-T = S imp =g(T/T K ), (1.32) 

where g is a scaling function which is universal for weak bare coupling. See Figure (^) . The behaviour 
of (7(2:) for small arguments is determined by RG-improved weak coupling perturbation theory. It's 
behaviour at low T is determined from the theory of the low energy fixed point. Its behaviour at 
arbitrary T/Tk is a property of the universal crossover between fixed points. It has been found from 
the Bethe ansatz. 

s 

imp 



In 2 - 




FIG. 8. Qualitative behaviour of the impurity entropy. 
Similarly, the susceptibility, at T — 0, is given by: 

X ' 



I b 
+ 



2tvvf Tk 

The ratio b/a is universal since the coupling constant (1/Tk) drops out. This is known as the Wilson 
Ratio. 

At high T, we must (for weak bare coupling) obtain approximately the results for a free spin: 



(1.33) 



X 



I 1 

+ ■ 



2-kvf 4T 

At lower T, using RG improved perturbation theory this becomes: 



/ 



X 



+ 



2ttv f 4T 



ln(T/lV 



+ . 



In general, we may write: 



^=x imp = ~/(T/:zv), 



(1.34) 



(1.35) 



(1.36) 



where /(T/Tjf) is another universal scaling function. See Figure m). 
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The temperature dependent part of the low T resistivity for the dilute random array is 2nd order 
in perturbation theory, 

P = Pu[l-d(^-) 2 ], (1.37) 

where d is another dimensionless constant. The second term comes from second order perturbation 
theory in the irrelevant coupling constant. Another universal ratio can be formed. I will discuss this 
in Sec. Ill, in some detail, using the CFT approach. 
We now expect a scaling behaviour: 

p(T) = nJ(T/T K ), 

sketched in Figure (px|). 




T 

FIG. 10. Qualitative behaviour of the resistivity. 



II. CONFORMAL FIELD THEORY TECHNIQUES 



It is very useful to bosonize free fermions to understand the Kondo effect. This allows separation 
of spin and charge degrees of freedom which greatly simplifies the problem. 
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We start by considering a left-moving spinless fermion field with Hamiltonian density: 



Define the current (=density) operator, 

J L {x - t) — : tp^tpL ■ (x,t) 

= ]im[ip L (x)ipL(x + e) - {0\^l(x)^ l (x + e)|0)] 

(Henceforth we generally drop the subscripts "L".) We will reformulate the theory in terms of 
currents (key to bosonization). Consider: 

J(x) J(x + e) as e — > 
= : tp^(x)tp(x) , ip^(x + e)ip(x + e) : 
+ [: i)\x)i>{x + e) : + : ^a:)^ + e) :]G(e) + G(e) 2 

G(e) = <0|^(ar)V + (a: + e)|0> = — . 

—it 

By Fermi statistics the 4-Fermi term vanishes as e — > 

: i/' t (a;)'i/)(x)V' t (a;)?/'(a;) : = — : t/^ (#)?/>''' (^^(a:)^^) : = 0. 
The second term becomes a derivative, 

lim[J(a;)J(a; + e) + ^] = lim — ^-[: ^(x)ip(x + e) : - : ^ (x + e) ip(x) :] 

= 2i : ip^ — ip : 
ax 

1 2 
TL = — J(x) + constant. 

Now consider the commutator, [J(x), J(y)]. The quartic and quadratic terms cancel. We must be 
careful about the divergent c- number part, 



[J(x),J(y)] = - 



d_ 

dx 



(x — y — iS) 2 (x — y + iS) 2 
1 1 



(5 0+) 



x — y — iS x — y + iS 



= 2ni — S(x — y) 
dx 

Now consider the free massless boson theory with Hamiltonian density (setting vf = 1): 



n 



i (dcj> 



2 \dt 



+ 



2 \dx) 



d 

[<t>{x), Qi^iv)] = iS ( x - v) 



We can again decompose it into the left and right-moving parts, 

(d t 2 - d x 2 )<f> = (d t + d x ){d t - d x )<p 

4>(x,t) = <j)h(x + t) + 4>r{x - t) 
(d t - d x )<t> L = d-4> L = 0, d+<t>R = 

h = \{d-^f + \{d + tf = \{d-<t> R f + \(d+<j> 

Consider the Hamiltonian density for a left-moving boson field: 

n = \{d + 4> L f 
[d+<t>L(x),d+cj>L(y)] = \4> + 4>'A + 4>'\ = 2 ^(* - y) 

Comparing to the Fermionic case, we see that: 
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since the commutation relations and Hamiltonian are the same. That means the operators are the 
same with appropriate boundary conditions. 

Let's compare the spectra. For the Fermionic case, choose boundary condition: 

^(0 = -V»(-0 (i.e. iPl(1) + iI>r{1)=01 fc=j(n+i), n = 0,±l,±2... (2.11) 

[See Figure (jE]). Note that we have shifted k by hp-] Consider the minimum energy state of charge 
Q (relative to the ground state). See Figure (|ll|). We have the single Fermion energy: 

E = v F k, (2.12) 



Q-l 



21 



(2.13) 



Now consider particle hole excitations relative to the Q-ground state: The most general particle-hole 
excitation is obtained by raising n m electrons by m levels, then n m _i electrons by m — 1 levels, etc. 
[See Figure @.] 



_ TVVF , 1 _ 2 \ 

E = -j-(t; Q +2^ n ™- m ) 



(2.14) 




k-k 



FIG. 11. A state with Q extra electrons added to the groundstate. 
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k-k 



FIG. 12. A particle-hole excitation in which three electrons are raised four levels and then one electron is raised three levels. 



Now consider the bosonic spectrum. What are the boundary conditions? Try the periodic one, 



</>(l) = 0H) => k : 



I 



(2.15) 



The m* single particle level has E The total energy is 



E 



TVVF 



(S^ n m ■ m), n„ 

1 



occupation number : 0, 1, 2, 



Where does the Q term in Eq. (2.14) come from? We need more general boundary condition on 



the boson field. Let be an angular variable: 

4>l(-1) = <Pl(1) + V^Q, <2 = 0,±1,±2, 



L(x + t) = ^l^--(x + t) + Y _i (e- iS r^ +t) a m + h.c), 

2 I, ^ ^/i I ^ m 



(2.16) 



(2.17) 



where a n 's are the annihilation operators and Q is the winding number, 



E : 



d.v 



2 Vat 



+ ■ 



2 \dx 



Here we have set vf = 1. We have the following correspondence: 

soliton <-> electron, oscillator <-> particle-hole pair. 

It is also possible to represent fermion operators in terms of the boson, 



ip L ~ e 



which gives the correct Green's function and implies the same angular definition of 4>l- 

For the Kondo effect we are also interested in the phase-shifted boundary condition: [See Figure 



(2.18) 



(2.19) 



iM0 = -hM-0. 



k — —n, (for fermions) 



E : 



TVVF 



0(0-1) 



E 



+ > n m m 



(2.20) 
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We have the degenerate ground state, Q = or 1, which correspond to an anti-periodic boundary 
condition on <f>, 

<t>{i) = 4>{-i) + V^(Q - \) 

E= i\ {Q -\ )2 + - = i { h Q{Q const - + (2 - 21) 

Now we include spin, i.e. we have 2-component electrons, 

H = iv F ip a] ^-ip a , (a = 1,2, summed). (2.22) 
ax 

Now we have charge and spin currents (or densities). We can write H in a manifestly SU(2) invariant 
way, quadratic in charge and spin currents: 

J=: V Qt ^ : , J = ^ a ^i3 (2.23) 



2 



Using: 



= 28% - 8%8] (2.24) 

J 2 — -- : ip ta ip a ^ ]l3 ipt3 ■ +7tV , " + -r~ipa +c-number, 
4 2 ax 

J 2 = : tf)^ a ijj a ijj^ tp3 : +2iip a+ — —ip a + c-number, 

ax 

H = ^-J 2 + ±- f 2 , (2.25) 
we have the following commutation relations, 

[J(x), J(y)] = 4:iri5'(x — y), (twice the result for the spinless case) 

[J{x), J z (y)] = \[Jr + Ji,Ji- Ji] = 0. (2.26) 

From [J, J] = 0, we see that H is sum of commuting charge and spin parts. 

[J a (x), J\y)\ = 2 7 n/> + [f °, |V ' *(* " ») + tr[§°> § - y) 

= 27vie abc J c {x) ■ 5(x -y)+ Tvi8 ab -^-8(x - y). (2.27) 

ax 

We obtain the Kac- Moody algebra of central charge k = 1. More generally the coefficient of the 
second term is multiplied by an integer k. 
Fourier transforming, 

J n = -L J dxe m f*J(x), [J%, J b m ] = ie abc J c n+m + \n8 ab 8 n ,- m (2.28) 

we have an co-dimensional generalization of the ordinary SU(2) Lie algebra. The spin part of the 
Hamiltonian is 



H s = yi ]T : J_„. J n : (2.29) 



The spectrum of H s is again determined by the algebra obeyed by the J n 's together with boundary 
conditions. The construction is similar to building representations of SU{2) from commutation 
relations, i.e. constructing raising operator, etc. 

In the k = 1 case we are considering here it is simplest to use: 

J{xf =3(J z (x)) 2 

= \((d + ^) 2 + (d + ^f) (2.30) 



15 



Now we have introduced two commuting charge and spin free massless bosons. SU(2) symmetry is 
now concealed but boundary condition on <j> s must respect it. Consider the spectrum of fermion 
theory with boundary condition: tj}(l) = —if>(—l), 



2 2 00 

2 + 2 + E m ^ + nL) 



2 2 

??l= — OO 



Change over to <f) c and <j> s . We can relabel occupation numbers, 



T I _ c s 



1 1 

-(x + t) + ... 



2^2 J 



, TT S Z , . 

+. = - T (, + t)+... 

Actually charge and spin bosons are not completely decoupled; we must require Q — 2S Z (mod 2), 
to correctly reproduce the free fermion spectrum. We see that the boundary conditions on <f> c and 
4> c are coupled. Now consider the phase-shifted case. 

7TVF ,1 f ^ .,s2 , rcz\2 



Redefine Q — 1 — > Q so 



E=^[ 1 -Q> + (Sr + ....], 



the same as before the phase shift, Eq. (2.32). One of the 0-energy single-particle states is filled, 
for Q = and there are 4 groundstates, 

(Q,5 z ) = (0,±i), (±1,0). 

Now Q = 2S Z + 1 (mod 2); i.e. we "glue" together charge and spin excitations in two different ways, 
either 

(even, integer) © (odd, half-integer) 
or (even, half-integer) © (odd, integer), 

depending on the boundary conditions. The -| phase shift simply reverses these "gluing conditions" . 

The set of all integer spin states form a "conformal tower" . They can be constructed from the Kac- 
Moody algebra by applying the raising operators J_ n to the lowest (singlet) state, with all spacings 
I2LE.. (integer). Likewise for all half-integer spin states, (s 2 ) 2 = 1+integer. Likewise for even and 
odd charge states. The K-M algebra determines uniquely conformal towers but boundary conditions 
determine which conformal towers occur in the spectrum and in which spin-charge combinations. 



III. CONFORMAL FIELD THEORY APPROACH TO THE KONDO EFFECT 



The chiral one-dimensional Hamiltonian density of Eq. QL19) and (1.14) is 



d a 
TL = —4>Lu + Ay4,™ -^-ipLp • S 5(x) (left-movers only) 

We rewrite it in terms of spin and charge currents only, 

H = -£-J 2 + ^{Jf + AJ • S 8{x). 

87T D7T 
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The Kondo interaction involves spin fields only, not charge fields: H = H 3 +H C . Henceforth we only 
consider the spin part. In Fourier transformed form, 



H s = y(g ^ ] J-n ■ Jn + A ^ ^ • S) 



r tq rb i ■ abc re I X a ^X 

[5 ,a ,5 ,i '] =i6 abc 5' c 
[S a ,4] = 



From calculating Green's functions for J(x) we could again reproduce perturbation theory d [* D = 

2 



-A 2 + • ••. That is a small A > grows. What is the infrared stable fixed point? Consider A — J 



where we may "complete the square" 



n — — oo 



[«^n "I" & i Jm + 5* ] — ie (>/,i+ m + S ) + — 5 5n,—m- 

H is quadratic in the new currents, J n = J„ + S , which obey the same Kac-Moody algebra! What 
is the spectrum of H(\ = |)? We must get back to Kac-Moody conformal towers for integer and 
half-integer spin. This follows from the KM algebra and the form of H (i.e. starting from the lowest 
state we produce the entire tower by applying the raising operators, J_„). 

Thus we find that the strong- coupling fixed point is the same as the weak- coupling fixed point. 
However, the total spin operator is now Jo = Jo + S. We consider impurity spin magnitude, s=l/2. 
Any integer-spin state becomes a 1/2-integer spin state and vice versa. 

Integer <-> 1/2-Integer. 

Presumably A = 2/3 is the strong coupling fixed point in this formulation of the problem, oo- 
coupling can become finite coupling under a redefinition, eg. 

^Lattice — ~, 3 



KM 



We expect the low-energy, large I spectrum to be KM conformal towers for any A. The effect of the 
Kondo interaction is to interchange the two conformal towers, Integer ^"Integer. See Figure 
(|l^). This is equivalent to a ^ phase-shift, 

(even, integer) © (odd,i-integer) 

t 

(even, ^-integer) © (odd, integer) 



Leading Irrelevant Operator: Specific Heat, Susceptibility, Wilson Ratio, Resistivity at T > 

At the stable fixed point S has disappeared; i.e. it is absorbed into J , 

J(x) = J(x) + 2nS S(x). 
What interactions could be generated in H e //? These only involve J ', not S. 

H s = ^-J(x) 2 + \ 1 J(0) 2 S(x). 

This is the only dimension-2 rotationally invariant operator in the spin sector. We have succeeded 
in reducing two dimension-2 operators to one. The other one is the charge-operator A2 J(0) 2 8(x), 
A2 = because there is no interaction in the charge sector (with other regularization we expect 
Ai ~ y— , A2 ~ -j, << Ai). 
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1/2 integer- s 



integer- s 



tower tower 

FIG. 13. At A = 2/3 the 1/2- integer- spin conformal tower is mapped into the integer-spin conformal tower. 

Now we calculate the specific heat and susceptibility to 1st order in Ai. 
Susceptibility of left-moving free fermions: 



0-th order A! = ^(n T - n{) = I J de v(e)[n(e + ^) - n(e - ^)] 



X = Y (for T « D) 



3T 

" 3T 2 ' 



■--^([J dxJ(x)] 2 J(0) 2 } + 



A simplifying trick is to replace: 



sn = \ij 2 (o)S(x) — ^ 2 M, 

which gives the same result to first order in A (only) by translational invariance of H at A = 0. Now 
the Hamiltonian density changes into 



We simply rescale H by a factor 



Equivalently in a thermal average, 



1 I 



x(Ai,T) = — <( / J)" > T(Al) 



l + 37rAl/Z 

n 37rAi 
[1-—]X0 

Z 3Ai 
2tt 2 ' 



x(0,T(Ai)) 
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where in the last equality the first term represents the bulk part and the second one, of order ~ , 
comes from the impurity part. Specific Heat: 

ttIT 

0-th order C = C C + C S , C c = C a = — . (3.16) 

Each free left-moving boson makes an identical contribution. 

d 

1st order in Ai C e (Xi, T) = — < H(Xi) > Al 

= C s (0,T(Ai)) 
ttI T 



3 1 + 3ttAi/Z 



nlT 



nXiT (3.17) 



The Wilson Ratio: 



SC, = _^ =2 SC* (3 . 18) 



Rw ^M*= 2= C (3.19) 



SC/C Cs 

measues the fraction of C coming from the spin degrees of free 

Doing more work, we can calculate the resistivity to 0(A 2 ).uO First we get the electron lifetime 
from the self-energy. The change in the 3D Green's function comes only from the ID s-wave part: 

G 3 (fi,f2) - G°(\fl - ra|) 

- ikFiri+r2) (G LR ( ri ,r 2 ) - G LR ,o(n,r2)) + h.c] 



8tt 2 



= G§(n)£G§(n,). (3.20) 

The self-energy E depends only on the frequency. It gets multiplied by the impurity concentation 
for a finite density (in the dilute limit). We must calculate the ID Green's function Glr{ti, r 2 , uj) 
perturbatively in A 

0(A?) : GLRin,^) = -Gl L (n,-r 2 ) 
= -G° LL (n + ra) 

= -Gl R (n,r 2 ), (3.21) 
where the (— ) sign comes from the change in boundary conditions, 

Gt R -G LR = -2Gl R + O(Xi) (3.22) 
To calculate to higher orders it is convenient to write the interaction as: 

J 2 = -f : V tc WV/3 : +|(V- ta ^^ - ^V>c) (3.23) 
To second order in Ai, we have the Feynman diagrams shown in Figure (|lj), giving: 

■£ R {u) = ^[2 + 37riAiW - i(37rAi)V - i(37rAi) 2 (u; 2 + tt 2 T 2 )]. (3.24) 
The first three terms give a phase-shift and the last term represents inelastic scattering. 
= ~2mJ^ ] + S inei (w) 

. TV 37rAlW 

5 = — | h ... 

2 2 



2/ m Eij(w) 

-[2- i(37rAi) V - i(37rAi) 2 (w 2 + tt 2 T 2 )] (3.25) 



r(o>) nv 2 4 
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The leading Ai dependence is O(Af) in this case. The O(Ai) term in E fl is real. We calculate the 
conductivity from the Kubo formula. (There is no contribution from the scattering vertex for pure 
s-wave scattering.) 



a(T) 



2e 2 


f d 3 k 


dn ' 


3m 3 J 


' (2tt)3 





k r(e k ) 



T(£fc) * 2^ [1 + I (37rAl )£ * + l (37rAl) (£ * + (7F T 



P(T) 



Srii 



-7T AiT ] 



All low-temperature properties are determined in terms of one unknown coupling constant Ai ^ 
Numerical or Bethe ansatz methods are needed to find the precise value of \i(D, A) oc ie 1 /*. 



(3.26) 



> x = ■ t 





FIG. 14. Feynman diagrams contributing to the electron self-energy up to second order in the leading irrelevant coupling 
constant of Eq. (3.8). 



IV. MULTI-CHANNEL KONDO EFFECT 



Normally there are several "channels" of electrons -e.g. different d-shell orbitals. A very simple 
and symmetric model is: 



//= E ^r^+^s- 



fe,Q!,i=l,2, . . . k 



k,k f a ,0% 



This model has SU(2) x SU(k) x (7(1) symmetry. Realistic systems do not have this full symmetry. 
To understand the potential ismplicability of this model we need to analyse the relevance of various 
types of symmetry breakingtil An interesting possible reicperimental application of the model was 
proposed by Ralph, Ludwig, von Delft and Buhrman.tJ In general, we let the impurity have an 
arbitrary spin, s, as well. 

Perturbation theory in A is similar to the result mentioned before: 



d\ 
d InD 

S 2 



- v )? + |^ 2 A 3 + 0[ks{s + 1)A 4 ] 
s(s+l). 



Does A — > oo as T — > 0? Let's suppose it does and check consistency. What is the groundstate for 
the lattice model of Eq. (1.2C), generalized to arbitrary k and s, at X/t — > oo? In the limit we just 
consider the single-site model: 



(4.1) 



(4.2) 



H = AS f -^|Vo, 



for A > (antiferromagnetic case) the minimum energy state has maximum spin for electrons at 
i.e. spin= k/2. Coupling this spin-k/2 to a spin-s, we don't get a singlet if s 7^ k/2, but rather 
an effective spin of size \s- k/2\. [See Figure ©.] The impurity is underscreened (k/2 < s) or 
overscreened (k/2 > s). 



(4.3) 
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FIG. 15. Formation of an effective spin at strong Kondo coupling, k = 3, s = 1 and s e g? = 1/2. 



Now let j « 1 be finite. Electrons on site ±1 can exchange an electron with 0. This gives an 
effective Kondo interaction: 

t 2 

A eff~X <<X 

See Figure (|l6|). What is the sign of A e g-? The coupling of the electron spins is antiferromagnetic: 
AgfrSei.o ■ Sei.i, with A e g- > (as in the Hubbard model). But we must combine spins 

S e ff = S + S e lfi- 

For | < s, S e ff \ | — S e i,o but, for | > s, S e f / 1 1 + S e i,o ■ So, ultimately, A e g- < in the underscreened 
case and A e g- > in the overscreened case. In the first (underscreened) case, the assumption A — > oo 
was consistent since a ferromagnetic A e g- — > under renormalizaton and this implies A — + oo, since 
A e fj ~ — j- 2 . In this case we expect the infrared fixed point to correspond to a decoupled spin of 
size s g g- — s ~ k/2 and free electrons with a 7r/2 phase shift. In the second (overscreened) case the 
oo-coupling fixed point is not consistent. Hence the fixed point occurs at intermediate coupling: 
This fixed point does not correspond to a simple boundary condition on electrons, instead it is a 
Non- Fermi-Liquid Fixed Point. See Figure ((lt|). 



| t I t 



FIG. 16. Effective Kondo interaction of 0(t 2 



o x 



c 

FIG. 17. RG flow of the Kondo coupling in the overscreened case. 



For the k=2, s=l/2 case we may think of the electrons (one from each channel) in the first layer 
around the impurity as aligning antiferromagnetically with the impurity. This overscreens it, leaving 
an effective s=l/2 impurity. The electrons in the next layer then overscreen this effective impurity, 
etc. At each stage we have an effective s=l/2 impurity. [See Figure (jl^).] Note in this special 
case that there is a duality between the weak and strong coupling unstable fixed points: they both 
contain an s=l/2 impurity. 
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FIG. 18. The overscreened case with s = 1/2, k — 2. 



A. Large-k Limit 



The /3-function is: 

/3 = A 2 -|a 3 + 0(A 4 ). 
If we only consider the first two terms, there is a fixed point at: 

A c w 2/fc. 

At this (small) value of A the quartic term, and all higher terms are 0(l/fc 4 ), whereas the quadratic 
and cubic terms are 0(l/k 2 ). Thus we may ignore all terms but the quadratic and cubic ones, for 
large k. The slope of the /3-function at the critical point is: 




This implies that the leading irrelevant coupling constant at the non-trivial (infrared) fixed point 
has dimension 2/k at large k, so that (A — A c ) scales as A 2 ^ fe . Thus the leading irrelevant operator 
has dimension (1+2/k). This is not an integer! This implies that this critical point is not a Fermi 
liquid. 



B. Current Algebra Approach 



We can gain some insight into the nature of the non-trivial critical point using the current algebra 
approach discussed in the previous section for the k=l case. It is now convenient to use a form 
of bosonization which separates spin, charge and flavour (i.e. channel) degrees of freedom. This 
representation is known as a conformal embedding. We introduce charge (J), spin (J) and flavour 
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(J A ) currents. A runs over the k 2 — 1 generators of SU(k). The corresponding elements of the 
algebra are written T A . These are traceless Hermitean matrices normalized so that: 

trT A T B = h AB , (4.8) 

and obeying the completeness relation: 



c <>a - T"a"c 



k 



A 

and the commutation relations: 



where the f ABC are the SU(k) structure constants. Thus the currents are: 

J = : i/) tlQ V»Q : 

7* / \ia i 

J = ip [ 



(4.9) 



[T A ,T B ]=if ABC T c , (4.10) 



r = r(n^. (4.ii) 



(All repeated indices are summed.) It can be seen using Eq. (4.9) that the free fermion Hamiltonian 
can be written in terms of these currents as: 

n = ik J2 + 2^2) P + Mk^Y) jAjA - ^ 

The currents J obey the SU{2) Kac-Moody algebra with central charge k and the currents J A obey 
the SU(k) Kac-Moody algebra with central charge 2: 

(4.13) 

The three types of currents commute with each of the other two types, as do the three parts of the 
Hamiltonian. Thus we have succeeded in expressing the Hamiltonian in terms of these three types 
of excitations: charge, spin and flavour. The Virasoro central charge c (proportietial to the specific 
heat) for a Hamiltonian quadratic in currents of a general group G at level k is:E3 

CG * - ZM^)Tfc' (4 - 14) 

where Dim (G) is the dimension of the group and Cv(G) is the quadratic Casimir in the fundamental 
representation. For SU(k) this has the value: 

C v (SU(k))=k. (4.15) 

Thus the total value of the central charge, c, is: 

C TOT = l + ^ + ^^ = 2,, (4.16) 

the correct value for 2k species of free fermions. Complicated "gluing conditions" must be imposed to 
correctly reproduce the free fermion spectra, with-yarious boundary conditions. These were worked 
out in general by Altshuler, Bauer and ItzyksonEJ The SU(2)k sector consists of k + 1 confp«pal 
towers, labelled by the spin of the lowest energy ("highest weight") state: s = 0, 1/2, 1, ...fc/2r 2 ll 33 l 
We may now treat the Kondo interaction much as in the single channel case. It only involves the 
spin sector which now becomes: 

ns = ^ T¥ / + XJ.§5( X ). (4.17) 

We see that we can always "complete the square" at a special value of A: 

where the Hamiltonian reduces to its free form after a shift of the current operators by S which 
preserves the KM algebra. We note that at large k this special value of A reduces to the one 
corresponding to the critical point: A c — > 2/k. 

While this observation is tantalizing, it leaves many open questions. We might expect that some 
rearranging of the (k + 1) SU(2)k conformal towers takes place at the critical point but precisely 
what is it? Does it correspond to some sort of boundary condition? If so what? How can we 
calculate thermodynamic quantities and Green's functions? To answer these questions we need to 
understand some more technical aspects of CFT in the presence of boundaries. 
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V. BOUNDARY CONFORMAL FIELD THEORY 



We will assume that the critical point corresponds to a conformally invariant boundary condition 
on the frae-jtheory. Using the general theory of conformally invariant boundary conditions developed 
by CardyEa we can completely solve for the critical properties of the model. Why assume that the 
critical point corresponds to such a boundary condition? It is convenient to work in the space- 
imaginary) time picture. The impurity then sits at the boundary, r = of the half-plane r > 
on which the Kondo effect is defined. If we consider calculating a two-point Green's function when 
both points are taken very far from the boundary (with their separation held fixed) then we expect 
to obtain bulk behaviour, unaffected by the boundary. [See Figure ([l9]).] This, at long distances 
and times is the conformally invariant behaviour of the free fermion system. Very close to the 
boundary, we certainly do not expect the behaviour to be scale invariant (let alone conformally 
invariant) because various microscopic scales become important. The longest of these scales is 
presumably the Kondo scale, £x ~ vf/Tl ~ ae 1 ^ uX . Beyond this distance, it is reasonable to expect 
scale-invariant behaviour. However, if the two points are far from each other compared to their 
distance from the boundary [Figure (po|)] then the behaviour is still influenced by the boundary 
even when both points are far from it. We have a sort of boundary-dependent termination of the 
bulk conformally invariant behaviour. The dependence on the details of the boundary (such as the 
value of £k) drops out. We may think of various types of boundaries as falling into universality 
classes, each corresponding to a type of conformally invariant behaviour. Rather remarkably, the 
above statements hold true whether we are dealing with a 2-dimensional classical statistical system 
with some boundary condition imposed, or dealing with a (l+l)-dimensional quantum system with 
some dynamical degrees of freedom living on the boundary. In fact, we already saw an example of 
this in the single-channel Kondo problem. The dynamical impurity drops out of the description of 
the low-energy physics and is replaced by a simple, scale- invariant boundary condition, tpL ~ —ipR- 




FIG. 19. The bulk limit. 




FIG. 20. The boundary limit. 
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Precisely what is meant by a conformally invariant boundary condition? Without boundaries, 
conformal transformations are analytic mappings of the complex plane: 



z = t + ix, 

into itself: 

z — ► w(z). 

(Henceforth, we set the Fermi velocity, vp — 1.) We may Taylor expand an arbitrary conformal 
transformation around the origin: 

oo 

w(z) = a n z n , 
o 

where the a n 's are arbitrary complex coefficients. They label the various generators of the conformal 
group. It is the fact that there is an infinite number of generators (i.e. coefficients) which makes 
conformal invariance so powerful in (1+1) dimensions. Now suppose that we have a boundary at 
x — 0, the real axis. At best, we might hope to have invariance under all transformations which 
leave the boundary fixed. This implies the condition: 

w(t)* = w(t). 



We see that there is still an infinite number of generators, corresponding to the a n 's of Eq. (5.3) 
except that now we must impose the conditions: 



We have reduced the (still oo) number of generators by a factor of 1/2. The fact that there is still an 
oo number of generators, even in the presence of a boundary, means that this boundary conformal 
symmetry remains extremely powerful. 

To exploit this symmetry, following Cardy, it is very convenient to consider a conformally invariant 
system defined on a cylinder of circumference f3 in the r-direction and length I in the x direction, 
with conformally invariant boundary conditions A and B at the two ends. [See Figure (|2l|).] From 
the quantum mechanical point of view, this corresponds to a finite temperature, T = 1//3. The 
partition function for this system is: 

Zab = tre _/3H ^, 

where we are careful to label the Hamiltonian by the boundary conditions as well as the length of 
the spatial interval, both of which help to determine the spectrum. Alternatively, we may make a 
modular transformation, r «-»■ x. Now the spatial interval, of length, f3, is periodic. We write the 
corresponding Hamiltonian as H p . The system propagates for a time interval I between initial and 
final states A and B. Thus we may equally well write: 

Zab =< A\e-' H p\B > . 



Equating these two expressions, Eq. (5.6) and (5/7) gives powerful constraints which allow us to 



determine the conformally invariant boundary conditions. 

h 1 




B 



FIG. 21. Cylinder of length I, circumference /3 with boundary conditions A and B at the two ends. 
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To proceed, we make a further weak assumption about the boundary conditions of interest. We 
assume that the momentum density operator, T — T vanishes at the boundary. This amounts to a 
type of unitarity condition. In the free fermion theory this becomes: 

V^Vl-MJ-V^'^M) =0. (5.8) 

Note that this is consistent with both boundary conditions that occured in the one-channel Kondo 
problem: ipL = ±tpR- 

Since T(t, x) = T(t + x) and fit, x) = f(t - x), it follows that 

f(t,x) =T(t,-x). (5.9) 

i.e. we may regard T as the analytic continuation of T to the negative axis. Thus, as in our 
previous discussion, instead of working with left and right movers on the half-line we may work with 
left-movers only on the entire line. Basically, the energy momentum density, T is unaware of the 
boundary condition. Hence, in calculating the spectrum of the system with boundary conditions A 
and B introduced above, we may regard the system as being defined periodically on a torus of length 
21 with left-movers only. The conformal towers of T are unaffected by the boundary conditions, A, 
B. However, which conformal towers occur does depend on these boundary conditions. We introduce 
the characters of the Virasoro algebra, for the various conformal towers: 

i 

where Ef(2l) are the energies in the aP 1 conformal tower for length 21. i.e.: 

Ef(2J) = js?-^, (5.11) 

where the xfs correspond to the (left) scaling dimensions of the operators in the theory and c is the 
conformal anomaly. The spectrum of H l AB can only consist of some combination of these conformal 
towers, i.e.: 

Zab =J2 n ABXa(e- 7 " 3/l ), (5.12) 

a 

where the n AB are some non-negative integers giving the multiplicity with which the various con- 
formal towers occur. Importantly, only these multiplicities depend on the boundary conditions, not 
the characters, which are a property of the bulk left-moving system. Thus, a specification of all 
possible multiplicities, n AB amounts to a specification of all possible boundary conditions A. The 
problem of specifying conformally invariant boundary conditions has been reduced to determining 
sets of integers, n a AB . For rational conformal field theories, where the number of conformal towers 
is finite, only a finite number of integers needs to be specified. 

Now let us focus on the boundary states, \A>. These must obey the operator condition: 

[T(x)-f(x)]\A>=0 (Var). (5.13) 

Fourier transforming with respect to x, this becomes: 

[L n - L n ]\A >= 0. (5.14) 

This implies that all boundary states, \A > must be linear combinations of the "Ishibashi states" @ 

| a >= [a; m > <g)ja; m >. (5.15) 



Here m labels all states in the a th conformal tower. The first and second factors in Eq. (5.15) refer 
to the left and right-moving sectors of the Hilbert Space. Thus we may write: 

\A>=^2\a>< aO\A> . (5.16) 

a 

Here, 



|o0 >= |o;0 > ®|a;0 >. (5.17) 
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(Note that while the states, \a; m > n > form a complete orthonormal set, the Ishibashi states, 
a > do not have finite norm.) Thus, specification of boundary states is reduced to determining the 
matrix elements, < aO\A >. (For rational conformal field theories, there is a finite number of such 
matrix elements.) Thus the partition function becomes: 

Zab = ^2 < A \ a0 >< a0 \ B >< a \ e ~' H " p \ a > • (5-18) 

a 

From the definition of the Ishibashi state, \a > we see that: 

<a\e- lH r\a>=Y, e ~ 2 ' KAm ' ( 5 - 19 ) 
m 

the factor of 2 in the exponent arising from the equal contribution to the energy from T and T. 
This can be written in terms of the characters: 

<a\e- lH p\a>= X a{e- M/f) ). (5.20) 
We are now in a position to equate these two expressions for Zab- 

Zab=J2 < A \ a0 >< a0 ' B > Xa(e- 47vl/l3 ) = n AB Xa(e-^ /l ). (5.21) 

a a 

This equation must be t cun f or all values of 1//3. It is very convenient to use the modular transfor- 
mation of the charactersE3£j 

b 

Here S b is known as the "modular S-matrix". (This name is rather unfortunate since this matrix 
has no connection with the scattering-matrix.) We thus obtain a set of equations relating the 
multiplicities, n a AB which determine the spectrum for a pair of boundary conditions and the matrix 
elements < aO\A > determining the boundary states: 

^2 S b n AB =< A\a0 >< aO\B > . (5.23) 

b 

We refer to these as Cardy's equations. They basically allow a determination of the boundary states 
and spectrum. 

How do we go about constructing boundary states and multiplicities which satisfy these equations? 
Generally, boundary states corresponding to trivial boundary conditions can be found by inspection, 
i.e., given n b AA we can find < a\A >. We can then generate new (sometimes non-trivial) boundary 
states by fusion, i.e. given any conformal tower, c, we can obtain a new boundary state \B > and 
new spectrum n a AB from the "fusion rule coefficients", N° b . These non-negative integers are defined 
by the operator product expansion (OPE) for (chiral) primary operators, (f> a . In general the (OPE) 
of (f> a with 0;, contains the operator <j> c N^ b times. In simple cases, such as occur in the Kondo 
problem, the N^ b 's a#er-all or I. In the case of SU(2)k, which will be relevant for the Kondo 
problem, the OPE isOE3 

j ® 3 = b'-/l. + !> \j -/I + 2, . . .,mm{j + j',k-j - ]}. (5.24) 

Note that this generalizes the ordinary angular momentum addition rules in a way which is consistent 
with the conformal tower structure of the theories (i.e. the fact that primaries only exist with 
j < k/2). Thus, 

= 1 (1 3 - i'| < j" < min{j + j' , k-j- j'} 

— otherwise. (5.25) 



The new boundary state, \B >, and multiplicities obtained by fusion with the conformal tower c 
are given by: 

nfl 

< aO\B > = < aO\A > -| 



an 

n°AB = ^2 N bc-n h AA- (5.26) 

b 
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Here labels the conformal tower of the identity operator. Importantly, the new boun dary state 



and multiplicities so obtained, obey Cardy's equation. The right-hand side of Eq. (5.23) becomes 



HQ 

< A\aO >< aO\B >=< A\aO >< aO\A > -|. (5. 
The left-hand side becomes: 

2J S b n AB = ^2 S b N dcn A A- (5- 
b b,d 



We now use a remarkable identity relating the modular S-matrix to the fusion rule coefficients, 
known as the Verlinde formulaO 

QCl QCL 



Y,SSNl = ^. (5. 



This gives: 



d 



proving that fusion does indeed give a new solution of Cardy's equations. The multiplicities, n% B 
are given by double fusion: 



^Nb c Na c n AA . (5. 



[Recall that \B > is obtained from \A > by fusion with the primary operator c] It can be checked 
that the Cardy equation with A = B is then obeyed. It is expected that, in general, we can generate 
a complete set of boundary states from an appropriate reference state by fusion with all possible 
conformal towers. 



VI. BOUNDARY CONFORMAL FIELD THEORY RESULTS ON THE MULTI-CHANNEL KONDO 

EFFECT 

A. Fusion and the Finite-Size Spectrum 

We are now in a position to bring to bear the full power of boundary conformal field theory on 
the Kondo problem. By the arguments at the beginning of Sec. V, we expect that the infrared 
fixed points describing the low-T properties of the Kondo Hamiltonian correspond to conformally 
invariant boundary conditions on free fermions. We might also expect that we could determine 
these boundary conditions and corresponding boundary states by fusion with appropriate operators 
beginning from some convenient, trivial, reference state. 

We actually already saw a simple example of this in Sec. Ill in the single channel, s — 1/2, Kondo 
problem. There we observed that the free fermion spectrum, with convenient boundary conditions 
could be written: 

(0, even) (1/2, odd). (£ 

Here and 1/2 label the SU (2)i KM conformal towers in the spin sector, while "even" and "odd" 
label the conformal towers in the charge sector. We argued that, after screening of the impurity 
spin, the infrared fixed point was described by free fermions with a 7r/2 phase shift, corresponding 
to a spectrum: 

(1/2, even) (0, odd). (f 
The change in the spectrum corresponds to the interchange of SU conformal towers: 

0^1/2. (e 

This indeed corresponds to fusion, with the spin- 1/2 prima ry fie ld of the WZW model. To see this 
note that the fusion rules for SU '(2)i are simply [from Eq. ( p.25| )]: 
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Thus for an s — 1/2 impurity, the infrared fixed point is given by fusion with the j = 1/2 primary. 
This is related to our completing the square argument. The new currents at the infrared fixed point, 
J , are related to the old ones, J, by: 

If J and J were ordinary spin operators, then the new spectrum would be given by the ordinary 
angular momentum addition rules. In the case at hand, where J and J are KM current operators, 
it is plausible that the spectrum is given by fusion with the spin-s representation, generalizing the 
ordinary angular momentum addition rules in a way which is consistent with the structure of the 



KM CFT. In particular, Eq. (3.5) implies, for half-integer s, that states of integer total spin are 
mapped into states of half-integer total spin, and vice versa, a property which follows from fusion 
with j = 1/2. 

This immediately suggests a way of determining the boundary condition for arbitrary number 
of channels, k and impurity spin magnitude, s: fusion with spin-s. Actually, while this is possible 
for s < k/2, corresponding to exact or overscreening, it is not possible in the underscreened case 
since there is no spin-s primary with which to fuse for s > k/2. Instead, in the underscreened 
case, we assume fusion with the maximal possible spin, namely k/2. This seems to correspond to 
the (in this case stable) strong coupling fixed point described in Sec. III. k/2 electrons partially 
screen the impurity. The fact that further screening is not possible is related to Fermi statistics. 
The maximal possible conduction electron spin state at the origin, for k channels is k/2. This is 
also essentially the reason why there are no primaries with larger spin, as can be seen from the 
corresponding bosonization of free fermions. We reiterate this essential point: The infrared fixed 
point in the k-channel spin-s Kondo problem is given by fusion with the spin-s primary for s < k/2 or 
with the spin k/2 primary for s > k/2. We have referred to this as the "fusion rules hypothesis". If 
the general assumption that the infrared fixed point should be described by a conformally invariant 
boundary condition is accepted, then this hypothesis starts to seem very plausible. The general 
method for generating new boundary conditions is by fusion. Since the Kondo interaction appears 
entirely in the spin sector of the theory we should expect that the fusion occurs in that sector. The 
current redefinition J — > J and various self-consistency checks all point towards this particular set 
of fusions. 

An immediate way of checking the fusion rule hypothesis, and more generally the applicability of 
the boundary CFT framework to this problem, is to work out in detail the finite size spectrum for 
a few values of k and s and compare with spectra obtained by numerical methods. 

Let us first consider the exactly screened and underscreened cases, s > k/2, where fusion occurs 
with the spin k/2 primary. In this case the fusion rules are simply: 

. ft k . 
J®2 = 2~ J - 

Each conformal tower is mapped into a uniquc-conformal tower. It can be shown that this gives the 
free fermion spectrum with a it/ 2 phase shift. u 

We demonstrate the case k=2, s=l in Tables (Q) and (|l|). Let us start with antiperiodic boundary 
conditions in the left-moving formalism: 

Ml) = -iM-0- 

Let us express this free fermion spectrum, for 2 spin components and 2 channels, in terms of 
products of conformal towers in the charge, spin and flavour sectors. In this case, the flavour sector 
corresponds to 577(2)2 as does the flavour sector. We will refer to the corresponding quantum 
numbers as j for ordinary spin and jf for flavour (or "pseudo-spin"). We need the energy of the 
"highest weight state" (i.e. groundstate) of each conformal tower. For the Kac-Moody conformal 
towers, therkighest weight state transforming under the representation R of the group G at level k 
has energy :c3 



E R n Cr 



I k + Ca ' 

where Cr is the quadratic Casimir in the R representation and A refers to the fundamental repre- 
sentation. For the case of SU(2) the representations are labelled by their spin, j and the Casimirs 
are: 
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Cj = j(j + 1). 



(6.9) 



We also need the energy for the charge sector. These can be worked out by generalizing the method 
used in the k=l case in Sec. II. The energy for the lowest charge Q excitation, for each species of 
fermion is: 

F _*Q 2 



as shown in Eq. (2.14). Altogether we obtain 4 terms like this for the 4 species of fermions (2 spin 



2 flavours). We can express the total energy in terms of the total charge: 

Q = Qn + Q12 + Q21 + Q22, 



and various difference variables. This gives: 



E 



191 
I 8 



+ . 



This gives the energy of the charge Q primary. In addition to the energy of the primary state we 
obtain additional terms in the energy corresponding to the excitation level in the charge, spin and 
flavour conformal towers: uq, n a and n/. These are non-negative integers. Altogether, we may 
write the energy of any state as: 



j'O' + i) , j/07 + 1) 



+ UQ + n s + rif 



For primary states, uq — n s = n/ = 0. Q must be integer and j and jf must be integer or half- 
integer. The allowed combinations of Q, j and jf are what we refer to as "gluing conditions". They 
depend on the boundary conditions. For antiperiodic boundary conditions, the allowed fermion 
momenta are 

k = 7r(ra + l/2)/J. 

The corresponding energy levels are drawn in Figure (|^). Note that the groundstate is unique. It 
has j = jf = Q = 0. Thus we must include the corresponding product of conformal towers in the 
spectrum. The single particle or single hole excitation has j — jf — 1/2 and Q = ±1. The energy 



(6.10) 



(6.11) 



(6.12) 



(6.13) 



(6.14) 



13/4 3/4 
8 4 4 



TV 1 

12' 



the right value. 2-particle excitations have lE/ir = 1, Q — 2 and either ha ve j = 1, jf — or j = 0, 
jf = 1, due to Fermi statistics. Again the energy is given correctly by Eq. (6.13). The lowest energy 



particle-hole excitations, with IE/it — 1, have Q = and various values of j and jf. In particular, 
they can have j — jf — 1. It turns out that these excitations are Kac-Moody primaries. Again 
the energy is given correctly by Eq. (6.13). It can be shownEj that these are all conformal towers 
that occur, except that we must allow arbitrary values of Q, mod 4. These conformal towers are 
summarized in Table (Dl). 



(6.15) 



TABLE I. Conformal towers appearing in the k = 2 free fermion spectrum with anti-periodic boundary conditions. 



Q (mod 4) 





±1 

2 
2 




1 

1/2 



1 



Jf 


1 

1/2 

1 







1 

1/2 

1 

1 
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Now consider fusion with the j = 1 primary. This has the effect of shuffling the spin conformal 
towers in the following way: 



-> 1 
1/2 -> 1/2 

1 0. (6.16) 

The spectrum of Table (|) goes into that of Table (Q) under this shuffling. It can be checked that 
this corresponds to free fermions with a 7r/2 phase shift, i.e. periodic boundary conditions, or a shift 
of the Fermi energy by 1/2 a level spacing, drawn in Figure Q). Now note that the groundstate 
is 2 4 = 16-fold degenerate, since the zero-energy level may be filled or empty for each species of 
fermion. The charge , Q, in Table (|fl|) is now measured relative to the symmetric case where 2 of 
these levels are filled and 2 are empty. Also note that if make the replacement: 

(6.17) 



we get back the previous spectrum of Table (p. Making this replacement in Eq. ( 3.13 ), we obtain 



El El Q 

7T 7T 2 

(ignoring a constant). This corresponds to shifting the Fermi energy by 1/2-spacing; i.e. a n/2 
phase shift. 



(6.18) 



TABLE II. Conformal towers appearing in the k = 2 free fermion spectrum with periodic boundary conditions, obtained 
from the Table I by fusion with j = 1. 



Q (mod 4) j jf (EI/tt) 

10 1/2 

11/2 

±1 1/2 1/2 1/2 

2 113/2 
2 1/2 
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In the overscreened case the fusion rules are more interesting. They lead to spectra which cannot 
be obtained by applying any simple linear boundary conditions to the free fermions. Thus we may 
refer to these as non-Fermi liquid fixed points. It might be possible to find some kind of non-linear 
description of the boundary conditions in this case. But note that a non-linear boundary condi- 
tion effectively introduces an interaction into the theory at the boundary. A boundary condition 
quadratic in the fermion fields might induce an additional condition quartic in fields, etc. Thus 
specification of non-linear boundary conditions could be very difficult. Cardy's formalism cleverly 
sidesteps this problem by the device of focussing on the boundary states instead of boundary condi- 
tions, and providing a method (fusion) for producing these boundary states. As was stated above, 
and we will continue to see in what follows, knowledge of the boundary states will determine all 
physical properties of the theory so nothing is lost by using this abstract description of the boundary 
condition. 

For the k = 2, s = 1/2 example, the fusion rules give: 

-> 1/2 
1/2 0© 1 

1 -> 1/2. (6.19) 
Now we get a larger number of conformal towers in the spectrum with this boundary condition, 



shown in Table (III). We have shifted the groundstate energy to 0, in this Table. Note that energies 
El /ir = 1/8 and 5/8 now occur. These do not correspond to any possible linear boundary conditions 
on free fermions. Note in particular that, due to particle-hole symmetry, only phase shifts of or 
7r/2 are allowed. These give half-integer energies, as we saw above and in Tables (|) and (|ll|). 

TABLE III. Conformal towers appearing in the k = 2 free fermion spectrum after fusion with j = 1/2. 

Q (mod 4) j j f (EI/tt) 

1/2 

1/2 1 1/2 

± 1 1/2 1/8 

± 1 1 1/2 5/8 

2 1/2 1/2 

2 1/2 11 
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This spectrum was compared with numerical work on the k = 2, s = 1/2 Kendo effect and 
the agreement was excellent (to within 5% for several of the lowest energy states) .EJ This provides 
evidence that the fusion rule hypothesis is correct in the overscreened case. 



B. Impurity Entropy 



We define the impurity entropy as: 

S imp (T) = limi^S^T) - So(l,T)], (6.20) 

where Sq(1,T) is the free fermion entropy, proportional to I, in the absence of the impurity. We 
will find an interesting, non-zero value for <Sj m p(0). Note that, for zero Kondo coupling, Sj m p = 
ln[s(s + 1)], simply reflecting the groundstate degeneracy of the free spin. In the case of exact 
screening, (k = 2s), <5j m p(0) = 0. For underscreening, 

Simp(0) = M*V + 1)], (6.21) 

where a' = s — k/2. What happens for overscreening? Surprisingly, we will obtain, in general, the 
log of a non-integer, implying a sort of "non-integer groundstate degeneracy" . 

To proceed, we show how to calculate Sj m p(0) from the boundary state. All calculations are 
done in the scaling limit, ignoring irrelevant operators, so that S- im p(T) is a constant, independent 
of T, and characterizing the particular bo unda ry condition. It is important, however, that we take 
the limit I —* oo first, as specified in Eq. ( |6.2C| ), at fixed, non-zero T. i.e. we are interested in the 
limit, 1/0 —> oo. Thus it is convenient to use the first expression for the partition function, Zab in 



Eq. (5.7) 



Zab = < A \ a0 >< a °\ B > Xa{e~ M/P ) -> e^ 76 " < A\00 >< 00\B > . (6.22) 



Here 1 00 > labels the groundstate in the conformal tower of the identity operator, c is the conformal 
anomaly. Thus the free energy is: 

Fab = -ttcT 2 1/6 - Tin < A\00 >< 00|B > . (6.23) 

The first term gives the specific heat: 

C = ttcTI/3 (6.24) 
and the second gives the impurity entropy: 

S imp = ln < A \°° >< °°\ B > ■ ( 6 - 25 ) 
This is a sum of contributions from the two boundaries, 

S- im p = Sa + Sb- (6.26) 

Thus we see that the "groundstate degeneracy" <?a, associated with boundary condition A is: 

exp[5 imp4 ] =< A|00 >= g A . (6.27) 

Here we have used our freedom to choose the phase of the boundary state so that <?a > 0. For 
our original, anti-periodic, boundary condition, g = 0. For the Kondo problem we expect the low 
T impurity entropy to be given by the value at the infrared fixed point. Since this is obtained by 



fusion with the spin-s (or k/2) operator, we obtain from Eq. (5.26) 



The modular S-matrix for SU(2)k is 



B0 



S=§- (6.28) 



^W^^sin 



7r(2j + l)(2j , + l) 



2 + fc 



(6.29) 
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/ i\ sin[7r(2 S + l)/(2 + fc)] 
g(s ' fc) = sin[^/(2 + fc)] ' (6 - 30) 

This formula agrees exactly with the Bethe ansatz result .0 This formula has various interesting 
properties. Recall that in the case of exact or underscreening (s > k/2) we must replace s by k/2 in 
this formula, in which case it reduces to 1. Thus the groundstate degeneracy is 1 for exact screening. 
For underscreening we must multiply g by (2s' + 1) to account for the decoupled, partially screened 
impurity. Note that, in the overscreened case, where s < k/2, we have: 

1 (6-31) 



2 + k 2 + k 2 + fc' 

so g > 1. In the case k — *• oo with s held fixed, g 2s + 1, i.e. the entropy of the impurity spin 
is hardly reduced at all by the Kondo interaction, corresponding to the fact that the critical point 
occurs at weak coupling. In general, for underscreening: 

Kff<2s + 1. (6.32) 

i.e. the free spin entropy is somewhat reduced, but not completely eliminated. Furthermore, g is 
not, in general, an integer. For instance, for k — 2 and s = 1/2, g — y/2. Thus we may say that 
there is a non-integer "groundstate degeneracy" . Note that in all cases the groundstate degeneracy is 
reduced under renormalization from the zero Kondo coupling fixed point to the infrared stable fixed 
point. This is a special case of what we believe to be a general result: the groundstate degeneracy, 
always decreases under renormalization. This appears to be related to Zamolodchikov's c-theoremEJ 
which states that the conformal anomaly parameter, c, always decreases under renormalization. 
The intuitive explanation of the c-theorem is that, as we probe lower energy scales, degrees of 
freedom which appeared approximately massless start to exhibit a mass. This freezes out their 
contribution to the specific heat, the slope of which can be taken as the definition of c. In the case 
of the "g-theorem" the intuitive explanation is that, as we probe lower energy scales, approximately 
degenerate levels of impurities exhibit small splittings, reducing the-idegeneracy. 

So far, only a perturbative proof of the g-theorem has been gii«n.li3 It is completely analogous to a 
perturbative proof of the c-theorem given by Cardy and LudwigJ2£l independently of Zamolodchikov's 
more general proof. For the g-theorem proof, we consider perturbing around a boundary CFT fixed 
point with a barely relevant boundary operator, i.e. the action is: 

S = S -\( dr<j>(0,T), (6.33) 
Jo 

where (f> has dimension 1 — y with < y « 1. The f3- function has the form: 

(3 = y\-bX 2 , (6.34) 

for some constant, b. There is a nearby fixed point at: 

A c = y/b. (6.35) 

It is possible to calculate the small change in g using renormalization group improved perturbation 
theory. This gives: 

Sg/g = -^V /36 2 < 0. (6.36) 



C. Boundary Green's Functions: Two-Point Functions, T=0 Resistivity 

In this sub-section we explain the basic concents for calculation of Green's functions in the presence 
of a conformally invariant boundary conditionO We then work out the case of two-point funtions 
in detail. Finally we show how this gives information about the Kondo problem. 

The most important point is the consequence of the identification of left and right-moving sectors, 
discussed in Sec. V. In general, in the bulk theory, a typical local operator is a product of left and 
right-moving factors: 

4>(x) = <j> L {x)4> R (x). (6.37) 
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Here x is the spatial co-ordinate; we suppress the time-dependence. However, in the presence of a 
boundary, we use: 

0h(x) = 4>l{-x). (6.38) 

Thus a local operator with left and right-moving factors becomes a bilocal operator with only left- 
moving factors: 

4>(x) -» 4> L {x)4> L {-x). (6.39) 

[See Figure (p2~|).] Thus a one-point function becomes a two-point function, two-point becomes 
four-point etc. 




FIG. 22. A local operator becomes effectively bilocal in the presence of a boundary. 



Henceforth, the number of points in the Green's function will refer to the larger number after 
this doubling due to the identification of left with right. In the remainder of this sub-section we 
show how to calculate boundary two-point functions. Four point functions are discussed in the next 
sub- section. 

Our bulk operators are normally defined so that < <f>(x) >bulk~ 0- ^or a semi-infinite plane 
with a boundary, this one-point function essentially becomes a 2-point function which may have a 
non-zero value: 

< 4>(x) >a=< H{x)4>l{-x) >a= j^sg. (6.40) 

Here d is the scaling dimension of (f>L, which does not depend on the boundary condition, A. On 
the otherphand, the coefficient, Ca does depend on the boundary condition. Following Cardy and 
Lewellen,LJ we may calculate Ca in terms of the boundary state \A >. We assume that is a 
primary field. 

This is done by making a conformal mapping from the semi-infinite cylinder to the semi-infinite 
plane: 

7TW 

z = itanh— . (6.41) 

P 

Writing: 

z = r + ix 

w =r + ix', (6.42) 

we see that, as x' goes from —(3/2 to /3/2, r goes from — oo to oo. [See Figure (p3|).] The semi-infinite 
plane is x > and the semi-infinite cylinder is t' > 0. Note that we are regarding x' (space) as 
the periodic variable on the cylinder, i.e. we are calculating imaginary-time propagation from the 
boundary state \A >. We will calculate the one-point function on the cylinder, first directly, then 
by obtaining it from the half-plane by conformal mapping. To obtain the correlation function on 
the infinite half-cylinder, it is convenient to take the limit of a finite cylinder, with boundary state 
00 > (the highest weight state of the identity conformal tower) at the other end. Thus, on the 
infinite cylinder, 

n , ,. <00| e -( T - T '» ff p^(0,0)e- T ' H p|^> faM . 
< </>(r ,0) >A=hm T ^ oa ■ -H- 2 ! • ( 6 - 43 ) 

< oo| e - TW £u > 
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Now we insert a complete set of states between <j> and \A >. Since is a primary field, 0|OO > 
gives only a sum of states in the conformal tower of <f>. For convenience, we also consider the limit, 
t' — > oo, so that only the highest weight state, |0O > survives. Thus, 



< 0(t ) > A ^ T ,^ C 



Note that 



< OO|0|0O >< <l>0\A > 
< 00\A > e 



E$ = (2n/P)2d. 



We need the matrix element < OO|0|0O >, for the periodic Hamiltonian Hp. This can be obtained 
from the Gpen's function < 00\<f>(r[ )4>{t2)\00 >, arising from a calculation on the infinite cylinder 
of radius /3£j This can be obtained by a conformal mapping from the infinite plane, giving: 



< OO|0(xO0(T2)|OO > 



P . , 7T , , 

-smh-(r 1 - t 2j 

TV p 



P ' 



< OO|0 00 > I e 



Thus we obtain: 



< OO|0|0 >= 



Thus the desired one-point function, on the half-cylinder is: 



< 4>{t') >, 



2tt 



< 4>^\ A > _(2x/5)2d-r' 

<00\A> £ 




23. Conformal mapping of a two-point function from the semi-infinite plane to the semi-infinite cylinde 



Now let us repeat this calculation, by conformal transformation from the semi-infinite plane. As 
argued above, the result on the half-plane takes the form: 



< 6{X) >A~- 



Ca 
(2x) 2 ' 



We now obtain the result on the half-cylinder by the conformal transformation of Eq. (3.41). This 
gives: 
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< T > 



1/2- cylinder 



dz . 2 d , , , plane 

7T 2 7rT ' 



c A 



. 7r sech 2 7rr'//3 



(2a:) M 

2d 



(3 2tanh ttt'//3 

2d 

2ttt'//3 



2tt 



Thus, equating Eq, ( 
denned in Eq. ( |s~4o| ) 



and (6.5C), we finally obtain the desired formula for the coefficient Ca, 
< 0O\A > 



Ca = 



< 00\A > 



If \A > is obtained by fusion with primary operator c from some reference state \F > (standing for 
free), then: 

r , ^free 

st/st 

~ A s°/s$ ' 



Eq. (3.40) can be immediately generalized to the case where the left and right factors occur at 
different points in the upper half-plane: 



< 4> L {z)4> R {z) > a = 



Ca 



[(t - r') + i(x + x') 



An important application of this formula in the Kondo problem is to the single fermion Green's 
function, < ipL{zyipR{z') >. In the case of periodic boundary conditions: 

< i>L a (x)i> RjP (x) > free = < i>l ia (x)i> LjP (-x) > 

2x ' 

We may obtain the one-point function at the Kondo fixed point by fusion. The fermion operator 
can be written as a product of spin flavour and charge operators: 

ipLja oc g a hje v 

Here g and h are the left moving factors of the primary fields of the WZW m odels , transforming 
under the fundamental representation, i.e. g has s=l/2. We may use Eq. (6.52). The Kondo 
boundary condition is obtained by fusion with the spin-s conformal tower, where s is the spin of the 
impurity. Thus: 



, i Ira ( \ i i \ . Kondo 



m si /2 /sy 2 



2x Ss j Sq 

cos[7r(2s + l)/(2 + fc)] ffig 
2x 



cos[tt/(2 + k)} 



Here we have used the SU(2) modular S-matrix given in Eq. (6.3C) . There are several interesting 
points to notice about this formula. First of all, we see th at for the case of exact or underscreening, 
where s = k/2, the cosine in the numerator in Eq. (6.56) becomes: 



cos7r(fc + l)/(fc + 2) = -cos7r/(fc + 2). 



Thus: 



. ; t / s. Kondo ^ it i 
< V>iWR > fc /2 = - < i> L i>R >f r , 

As expected, this corresponds to a n/2 phase shift: 

MO) = e 2tS i> L (0) = -^l(O). 
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In general, we may define a one-particle into one-particle S-matrix element , by: 



< 1p\lpR >= S 



;1) < >free> 



with: 



(!) _ COs[7T(2s + l)/(2 + fc)] 



gw = 



)s[tt/(2 + k)] 



We see that in the overscreened case, IS 1 *- 1 ^] < 1. S (1 ^ is the matrix element, at the Fermi energy, 
for a single electron to scatter off the impurity into a single electron. In the case where IS 1 ' 1 '! = 1, 
we see, by unitarity, that there is zero probability for a single electron to scatter into anything but 
a single electron at the Fermi energy, at zero temperature. This is precisely the starting point for 
Landau's Fermi liquid theory. Thus we refer to such cases as Fermi liquid boundary conditions. In 
the overscreened case, where |5 | < 1, this inelastic scattering probability is non-zero so we have 
non-Fermi liquid boundary conditions. It is interesting to note that in the large k limit (with s held 
fixed), S 1 ^ — > 1, corresponding to the fixed point occuring at weak coupling. We also note that, 
for k = 2, s = 1/2, = 0. This is related to the symmetry between the zero coupling = 1) 

and infinite coupling (S^ = —1) fixed point, mentioned above. In this case the one particle to one 
particle scattering rate vanishes! 

From the electron self-energy we caru-Qbtain the lifetime and hence the resistivity for a dilute 
array of impurities by the Kubo formulali£l giving a T = resistivity: 



P(0) 



3ni 



kit(evvF) 2 



S 



(i) 



Here rij is the impurity density. The first factor is the "unitary limit", i.e. this is the largest possible 
resistivity that can occur for a dilute array of non-magnetic impurities. It is only realised in the 
Kondo problem in the Fermi liquid case (exact or overscreened). Otherwise the resistivity is reduced 



by the factor 



which goes to zero for large k. 



D. Four Point Boundary Green's Functions, Spin-Density Green's Function 

In this sub-section we sketch how four-point functions of chiral operators (or two-point functions 
of non-chiral ones) are calculated in the presence of a boundary, using the particular example of 
the fermion four-point function with the Kondo boundary condition. For more details, see Ref. ( 
[l6| ). Similarly to the case of the two-point function, we will find that the four-point function is 
determined by bulk properties up to a constant which can be expressed in terms of matrix elements 
involving the boundary state. 

We consider the Green's function: 

G=< ^L a i(zi)i> t L J {z2)i}Rf3j{z 3 )i}li'''(z 4 ) > . 

We suppress spin and flavour indices in what follows, but it must be understood that we are dealing 
with tensors throughout. The first step is to regard the right-moving fields as reflected left-moving 
ones: 

tJjr(z) — ► ip L (z*). 

to express G as a product of charge, spin and flavour Green's functions: 
G = GcGsGf. 

The form of G c is unique up to a multiplicative constant, and is unaffected by the boundary. G s 
and Gf are partially determined by the bulk conforrpal field theories, i.e. they must be solutions of 
the linear Knizhnik-Zamolodchikov (KZ) equations.c3 These equations have two solutions in both 
cases. We label them G s (p) and Gf(q) where p and q take values and 1. They are tensors in spin 
and flavour space; we suppress these indices. Thus we may schematically write the solution as : 

i 

G =J2 ^P,iGcG 3 (p)Gf{q). 

p,q = 



Then we use Eq. ( 3.55 



38 



The four constants, a p , q remain to be determined. It turns out that three of them can be determined 
from general considerations, with only one (oi,i) depending on the boundary conditions. 

To see this, and to understand better how the solutions G s {p) and Gflq) are defined, it is conve- 
nient to consider the limit Z\ z 2 and 23—5-24, illustrated in Figure (b4|). In the limit we may use 
the operator product expansion of tl>i,(z\) with ifr L {z^). Since these points are at a fixed distance 
from the boundary as they approach each other, the bulk OPE applies. This is simply the trivial 
OPE of free fermions: 

i>( Zl )^{z 2 ) — + J + J+J A + O a s d Of + ... (6.67) 

zi - 22 

Here we have used the fact that the spin and flavour currents are bilinear in the fermion fields. 
and 0°j d are the primary fields in the adjoint representation of the spin and flavour groups. These 
have scaling dimension 2/(2 + k) and k/(2 + k) respectively, so their product has scaling dimension 
1 and corresponds to the product of fermion fields with spin and flavour traces subtracted. Note, 
importantly, that neither adjoint primary can appear by itself, since it has a fractional scaling 
dimension that does not occur in the free fermion OPE. 

• z 

z* 2 

1 



z* • • z* 



FIG. 24. The bulk limit, zi -» z 2 , 23 -> 24. 



Now let us consider the bosonized expression for the fermion fields, of Eq. ( p.55[ ). Consider the 
OPE of each factor separately: 

e I ^ 7I ^(2 1 )e- l ^ 7I ^(2 2 ) ~ (2! - Z 2 )' 1/2k + ... 

g{zi)g\z 2 ) ~ (zi - Z2 )- 3 ^ 2+k » + (zi - ^/W+^of + ... 

h{zi)h\ Z2 ) ~ (Zl - 22 )-(* 2 -D/[M2+k)] + {zi _ Z2) l/M2+k) ]o a d + (g 6g) 

A general solution of the KZ equations for the four-point spin and flavour Green's functions, G 3 and 
Gf will have singularities corresponding to both the singlet and adjoint terms in the OPE shown 
above. The two independent solutions, G s (p) and Gf(q) can be chosen so that only the singlet 
appears for p = and q = and only the adjoint appears for p — 1 and q = 1. It is now clear that, 
in order for G to reproduce free fermio n beh aviour in the limit 21 — *■ z 2 , we must demand that the 



coefficients ai,o and 00,1 vanish in Eq. ( 6.66 ) so that the spin or flavour adjoint field doesn't occur 
without being multiplied by the other adjoint field. Similarly, requiring the correct normalization 
for the singlet- singlet singularity in this limit: 

G^-, — 1 w . 1 , v (6-69) 

(21 -22) [zl - zl) 

determines the non-zero value of the coefficient ao,o- i.e. it has the same value as in the free fermion 
case, independent of the boundary conditions. Thus only ai,i remains to be determined by detailed 
consideration of the boundary condition. 
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From considering the limit zi — > z-z, 23 — » 24, we see that: 

01,1 oc< [of of ](*i)[o; d of >a . (6.70) 

This is a two-point Green's function with the two points staddling the boundary. Thus its normal- 
ization does depend on the particular boundary condition, A. In fact this n orma lization is precisely 



the coefficient Ca which was c alcuated in the previous sub-section, Eq. (6.52). Thus we obtain 



upon identifying <f> in Eq. (3.52) with the adjoint (spin I) primary in the spin sector, for the Kondo 
boundary condition: 



^s/^o „rree i a V1 \ 



Si/ Sp free 
S2/S* 

This result completes the determination of the four point function. The general form of G is rather 
complicated, involving hypergeometric functions. We only give explicit results here for the simplest 
non-trivial case, k = 2, s — 1/2. We consider only the most singular (2ttF part) of the spin density 
Green's function. Writing: 

S a {f,r) = e 2 ^ r ^ L ^ R + h.c, (6.72) 

this becomes: 

1/2 



<S a (ri,T 1 )S b (f 2 ,T 2 )>- ' '' 



2cos2A:F(ri + ri) + [ 2 + - ] cos2fcF(ri - r 2 ) 



FIG. 25. The boundary limit z\ — * zf, z<i — * zf. 

This change at the boundary can be understood as follows. As explained in the previous sub- 
section, in the presence of a boundary, the spin density becomes a bilocal operator, at the points z 
and z* . As r — *■ we use the O.P.E. to express it as a sum of local operators. In general, any operator 
in the O.P.E. of the two left-moving operators Ol and Ol may appear. The coefficients of the various 
terms in the O.P.E. (some of which may be zero) depend on the particular boundary condition. This 
is very natural since we are taking an O.P.E. for two points which straddle the boundary. [See Figure 
(0).] In this case we are obtaining the spin adjoint primary field in the O.P.E., unaccompanied by 
the flavour adjoint, as it is in the bulk O.P.E. We give a general prescription for determining the 
boundary operator content in the next sub-section. 



(6.73) 



8n 4 rjr?, \zi - z 2 *| 2 

Here the ... represents less singular terms and r\ is the cross-ratio: 

(*i - - zj) _ ^2 , , 

1- (zi - z$)(z 2 - z{) (n - r 2 ) 2 + (n + r 2 ) 2 ' k ' ' 

Note that in the bulk limit, \z\ — z%\ « r\,r%, [Figure ( pif ) with z\ = Z3 and Z2 = Z4] 77 — > 1 and 
the Green's function reduces to its bulk form: 

. 5,^ , 8 ab cos2fc F (ri - r 2 ) 

< S {ri,n)S (r 2 ,T 2 >^ 0422 7 \T~n \2' ( 6 - 75 ) 

87r 4 r 2 r2 (n — r 2 ) 2 + (n — r 2 ) 2 

the free fermion result. On the other hand, in the boundary limit, [Figure (pj|)1 T\,ri « \ti — t 2 \, 
r] — * and 

. n c 6/^ <5 ab cos2fcFricos2fc F r2 , . 

< S (ri,n)5 (r 2 ,r 2 >^ —— -. (6.76) 

47r 4 (rir2)' i/2 |ri — T2 

The r-dependence implies that S a has scaling dimension 1/2 at the boundary, although it has 
dimension 1 in the bulk, corresponding to the spin current operator. 
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E. Boundary Operator Content and Leading Irrelevant Operator: Specific Heat, Susceptibility, Wilson Ratio, 

Resistivity at T > 



As we saw in Sees. I and III, the leading irrelevant operator plays a very important role in the 
Kondo problem, determining the low temperature behaviour of the specific heat, susceptibility and 
resistivity. One of the novel features of a non-Fermi liquid boundary condition is that boundary 
operators may appear which do not occur in the bulk theory. As explained in the previous sub- 
sections, this is a consequence of the bilocal nature of operators in the presence of a boundary and the 
fact that the O.P.E. coefficients depend on the particular boundary condition. In this sub-section we 
derive a general formula which gives all boundary operators that occur with a particular boundary 
condition. Then we analyse the effect of the leading irrelevant operator in the overscreened Kondo 
problem. 

We can identify the boundary operator content from a general relationship between the finite-size 
spectrum and operator content. This is established by a conformal mapping from the semi-infinite 
plane to the infinite strip. We consider a correlation function for some primary operator, O, on the 
semi-infinite plane: 

<0(r 1 )O t (r 2 ) > A ■ 



[See Figure (|2(j).] Now we make the conformal mapping: 

j 7TW II 
V in < 



Here w is on the strip: 



We define: 



-1/2 < lmw < 1/2. 



z = r + ix 
w = u + iv. 



It is convenient to assume ti,ts > so that Vi = v 2 = 0, as shown in Figure ([26|). The correlation 
function on the infinite plane has the form: 

< 0(71)0* fa) >A~- 



(Tl — T 2 ) 2x ' 

From this we obtain the correlation function on the strip: 

H(«0f§(«2) 



<0{ Ul )0\u 2 ) >AA = 



[z(lll) - z(u 2 )} 2 

— sinn —(wi - m 2 ) 

7T Zt 



Here AA denotes the correlation function on the strip with boundary condition A on both sides. As 
U2 — mi — > <x this approaches 



<0(r,)O t (r 2 )>^ 



-ttx(u2 — ui)/1 



Alternatively, we may evaluate the correlation function on the strip by inserting a complete set of 
states: 

< 0{ Ul )0\u 2 ) > AA = I <0\O\n > AA |2 e - B " c " 2 - ul) . 



(6.77) 
(6.78) 
(6.79) 

(6.80) 

(6.81) 



(6.82) 



(6.83) 



(6.84) 



Here we get the eigenstates on the strip with boundary conditions of type 'A' on both sides. As 
M2 — Mi — > oo, the lowest energy state created from the groundstate by O dominates. This is simply 
the highest weight state corresponding to the primary field O. Clearly, the corresponding energy is: 

E n = nx/l. 

Thus we see that the primary boundary operators with boundary condition A, are in one-to corre- 
spondance with the conformal towers appearing in the spectrum with boundary conditions AA. 



(6.85) 
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FIG. 26. Conformal mapping of a boundary two-point function from the semi-infinite plane to the finite strip. 



Actually we have made an assumption here, that the operator O does not change the boundary 
condition. This appears to be reasonable for the Kondo problem, but in some situations, such as 
the X-ray -edge singularity, it is neccessary to consider more general boundary condition changing 
operators.LJ 

Thus to obtain the boundary operator content with a boundary condition A, we simply need 
to know the finite size spectrum with the boundary condition A at both sides. In the case of the 
Kondo boundary condition, this is obtained by "double fusion", i.e. beginning with the free fermion 
spectrum, fusing once with the spin-s primary gives the spectrum with free fermion boundary 
conditions at one end and Kondo boundary conditions at the other. Fusing a second time with the 
spin-s primary gives the spectrum with Kondo boundary conditions at both ends, as mentioned in 
Sec. III. The fact that double fusion gives the spectrum with Kondo boundary conditions at both 
ends can also be seen by a "completing the square" argument, i.e., if we have Kondo impurities at 
x = and x — I, then at a special value of both Kondo couplings, Ai = A2 = 2/3, the Hamiltonian 
returns to its free form if we redefine the SU (2) currents by: 

J n = J„ + S ? i + (-l) n 5 f 2 . (6.86) 

Double fusion represents a generalization of ordinary angular momentum addition rules, for the case 
of two impurity spins, which is consistent with the conformal tower structure of the WZW model. 
With free fermion boundary conditions, the boundary operators are simply the bulk free fermion 
operators. After double fusion, new boundary operators, not present in the free fermion theory, 
are generally present. However, an exception occurs in the case s = k/2 corresonding to exact or 
underscreening. In this case under double fusion the jth conformal tower is mapped into itself: 

3 -> k/2 - j -» j. (6.87) 

Thus the free fermion operator content is recovered, corresponding to a Fermi liquid boundary 
condition. 

On the other hand, in the overscreened case, k > 2s, non-trivial operators are always obtained. 
In particular, the spin j = 1, charge Q = flavour singlet operator, <f>, always occurs. This follows 
from double fusion with the identity operator: 

(Q = 0, j = 0, flavour singlet) — > (Q = 0,j = s, flavour singlet) 

(Q = 0,3 = 0, 1, ...min(2s,fc - 2s), flavour singlet). (6.88) 

For overscreening, min(2s, k — 2s) > 1, so <f> (with j = I) always occurs. <j> itself, cannot appear in 
the effective Hamiltonian, since it is not s spin singlet. However, its first descendent, J_i • tf> can 
appear since it is invariant under all the symmetries of the underlying Kondo Hamiltonian. It has 
dimension, 

x = 1+ 2Tk- (6 ' 89) 

This dimension is > 1 meaning that the operator is irrelevant. On the other hand, the dimension is 
< 2 so it is more relevant than J 2 , the leading irrelevant operator inpfche Fermi liquid case. It can 
be seen to be the leading irrelevant operator in all overscreened casesJiil 
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We now discuss the calculation of various quantities to lowest non- vanishing order in the leading 
irrelevant term in the effective Hamiltonian: 

SH = A1J-1 • $. 

We begin with the specific heat. The simple change of variables trick that we used in Sec. Ill for the 
Fermi liquid case doesn't work here. In fact, since J_i • <j> is a Virasoro primary operator (although 
a Kac-Moody descendant) its one-point function continues to vanish even at finite temperature. 
Consequently, the leading contribution to the specific heat is second order in Ai. The contribution 
to the impurity free energy of 0(X{) is: 



-PA 



imp 



0/2 



0/2 



0/2 



dr 2 T < J_i ■ 4>(ti)J-i ■ <j>(T 2 ) > 



0/2 



Since J_i ■ <j> is a Virasoro primary, this two-point function is given by: 

3(2 + fc/2) 



< J- 



(n) J_i • <P(t 2 ) > 



Ti - r 2 



2(1+A) 



3(2 + fc/2) 



Ifsin^n-^I^+A)' 



Here 



A = 



2 + fc' 

and we have assumed that d> has unit normalization: 



< cb a (nH b (r 2 ) >-- 



\Tl — T 2 



The normalization of J_i • <f> can then be fixed using the Kac-Moody algebra. A naive rescaling 
argument implies that: 



/imp ^ 



1+2A 



Actually the integral in Eq. (6.91) requires an ultraviolet cut-off. However, this only leads to 
additional terms in /j m p that vanish more rapidly as T — *■ 0. Evaluating the universal term in Eq. 
(6.91) explicitly, we obtain: 



d 2 h 



C- — —T 
°imp A q T 2 

(for k > 2) 



imp _ 2 



T^7T 1+ ^(2A^3(-+2)i 



?or k , o^ 1^ ( 1 /2- A)r(l/2) 



r(i-A) 



+ 0(T) 



+ o(\It 3A ) + .. 



Here T(x) is Euler's Gamma function. Note that this result is independent of the magnitude of the 
impurity spin, s. Also note that since, for k > 2, A = 2/(2 + k) < 1/2, so the impurity specific heat 
is more singular than in the Fermi liquid case. In the particular case k = 2, where A = 1/2, we 
obtain from Eq. (6.91): 

Cimp = A?7r 2 9Tm(T A -/T) + ... 

Note that this is only more singular by a log than in the Fermi liquid case. 
The impurity susceptibility is also second order in Ai, 



Ximp oc / dT 1 dT 2 dx 1 dx2 < J-i ■ <?!>(ti)J_i • 4>(r 2 )J (Q,xi)J (0,x 2 ) > 
Again, by a scaling argument we see that: 

Ximp x A i J 



By evaluating the integral in Eq. ( B.9S ) explicitly, we can caluculate the Wilson ratio, Rw for the 
overscreened Kondo problem: 
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XimpC bulk _ (2 + fc/2)(2 + fc) 2 

tiw ~ 



C'imp Xbulk 18 

We can calculate the Wilson ratio in the exactly screened, Fermi liquid, case, much more simply. 
In this case, the simple rescaling argument of Sec. Ill always works giving: 

Cimp/^ _ 



^imp/X 

Here C s is the part of the bulk specific heat coming from the spin degrees of freedom. Hence the 
Wilson ratio is given by 

It simply measures the ratio of the total bulk specific heat to that coming from the spin degrees of 
freedom. Using: 

where c is the conformal anomaly, and the fact that c — 2k for 2k species of free fermions and 
c = 3k/ (2 + k) for the level k SU(2) WZW model, we obtain: 

2(2 + fc) 
= jj , 

for the exactly screened Kondo problem. 

Finally we consider the T-dependence of the impurity resistivity. For more details, see Ref. ( [j^) . 
We may write the electron Green's function in the one-dimensional theory with the impurity at the 
origin in the form: 

G(n, — ra,w n ) — Go(n + r 2 ,Un) = Go(n, w„)£i(w„)Go(»"2, w„). 

In the three dimensional theory with a dilute random array of Kondo impurities, of density m the 
self-energy becomes: 

v 



The electron life-time is given by: 



- = -2ImZ£, 

T 



where the superscript R denotes the retarded Green's function, obtained by analytic continuation 
from the Matsubara Green's function. The conductivity is then expressed in terms of the lifetime 
in the standard way: 



2e 2 f d 3 p 

a — 



3m? J (2tt) 3 



-dni 



P r( £ p) 



Thus, we just need to calculate the one-dimensional self-energy, Si. This was already done, at T=0, 
in Subsection (VIC). There we showed that: 

G = G S W , 

wher e is the scattering matrix element for 1-electron to go into 1-electron, given explicitly in 
Eq. (6.61). This gave the resistivity: 

1 - S« 

P=^— PU, 

where pu is the unitary limit resistivity. To obtain the leading T-dependence, we do perturbation 
theory in the leading irrelevant coupling constant, Ai- In this case, there is a non-zero first order 
contribution: 



<5Ei oc Ai < ipJ-i ■ W) 1 > 
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Since Ai has scaling dimension A, it follows that: 

1 - 

p = pu [l+aAiT A + ...] , (6.112) 

where a is a dimensionless constant which can be obtained by explicit evaluation of first order 
perturbation theory. After a rather long calculation, a can be expressed in terms of an integral 
over hypergeometric functions. We evaluated this integral numerically for the case k = 2, where 
A = 1/2, obtaining: 

a = 4v^F. (6.113) 

Thus we can form another universal ratio involving the square of the temperature-dependent term 
in the resistivity, and the specific heat coefficient. 

Note that in the Fermi liquid case the resistivity had the form: 

p = pu[l - a\(T 2 ]. (6.114) 

The leading T-dependence is second order in Ai and T. In the non-Fermi liquid case the sign of the 
leading T-dependence depends on the sign of the irrelevant coupling constant, Ai. As we increase the 
Kondo coupling constant, so as to pass through the non-trivial critical point, at Kondo couplings of 
O(l), the leading irrelevant coupling constant, Ai, should change sign. Thus the sign of the leading 
T-dependent term in the resistivity actually switches at the critical point. We can't determine the 
value of Ai from our methods. However, if we make the plausible assumption that Ai < for weak 
Kondo coupling and hence Ai > for strong Kondo coupling, then the resistivity behaves, at low 
T as shown in Figure (^) (since the constant, a > 0). For weak Kondo coupling the resistivity 
decreases with T, exhibiting a power law singularity at T = 0. On the other hand, for strong Kondo 
coupling, the resistivity initially increases with T [Figure J27|)]. This is quite reasonable for very 
strong Kondo coupling where, at high T the effective renormalized Kondo coupling is close to the 
unstable fixed point at oo. Hence, we expect to obtain the unitary limit resistivity at high T in this 
case. This requires p to increase with T. 




« 1 c (1) 
P 1-S — 

u 



large X 



FIG. 27. Qualitative behaviour of the resistivity in the overscreened case. 



This vT behaviour of the resistivity, and other scaling behaviour predicted by the present ap- 
proach for the case k — 2, s — 1/2, were observed in recent experiments. E3 However, the intep«ifet|e 
tion of these experiments in terms of the 2-channel Kondo problem is controversial at present J13'cd'E 
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